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Abstract 

Brooding over pions, wave-packets and Bose-Einstein correlations, we present a recently 
obtained analitical solution to a pion laser model, which may describe the final state of pions 
in high energy heavy ion collisions. 



' Introduction. At the planned RHIC accelerator pions, the lightest observable strongly interact- 

Q ■ ing particles, will be produced in thousands in a unit interval of the longitudinal phase-space. At 

^ . such high multiplicities, the bosonic nature of the pions and the resulting Bose-Einstein statistics 

! may dominate the behaviour of the observable physical quantities, like the multiplicity distribu- 

G^ I tion, the single-particle momentum spectra and the second order intensity correlation function 



(frequently referred to as the Bose-Einstein correlation function or BECF). Present experiments 
at CERN SPS involving central collisions of two almost completely stripped nuclei (heavy ions) 
^ . already produce cca 600 similar pions in a unit value of rapidity, when Pb + Pb reactions are 
^ ! measured at 160 AGeV laboratory bombarding energy. At the RHIC collider, to be completed by 
qh| 1999, thousands of pions could be produced in a unit rapidity interval in a 40 TeV collision of Au 
nuclei |jl|, ^ . This multiplicity is large enough for the higher order Bose-Einstein symmetrizations 
to reveal their effects on the observables. We report here about a simple model |0, Q that is ana- 
lytically solvable both in the dense and the rare Bose-gas limit. The model is based on stimulated 
emission of pion wave-packets that is induced by the overlap of the multi-pion wave-packet states. 
The presentation closely follows the lines of ref. |^. 

Specification of the Pion-Laser Model. If the number of pions in a unit value of phase-space is 
large enough these bosons may condense into the same quantum state and a pion laser could be 
created Similarly to this process, when a large number of bosonic atoms are collected in a 
magnetic trap and cooled down to increase their density in phase-space, the bosonic nature of the 
atoms reveals itself in the formation of a Bose-Einstein condensate 0, a macroscopic quantum 
state. Such a condensation mechanism may provide the key to the formation of atomic lasers in 
condensed matter physics and to the formation of pion lasers in high energy particle and heavy 
ion physics, reviewed recently in refs J^, |^. 

The density matrix of a generic quantum mechanical system is 

oo 

P=I]P"/3„, (1) 

n=0 



Here the index n characterizes sub-systems with particle number fixed to ra, the multiphcity 
distribution is prescribed by the set of {Pn}^o' normahzed as I]^oP" ~ ^- '^^^ density matrixes 
are normahzed as Trp = 1 and Trp„ = 1, where 

Pn=jdai...dan Pn(ai, •••,«„) |ai, ...,a„)(ai, (2) 

and the states |ai, ...,«„) denote properly normalized n-particle wave-packet boson states. 
A wave packet creation operator is 



«I = «nP), (3) 



TTcr^ 



where = (^j, vTj, (jj, tj) refers to the parameters of the wave packet i: the center in space, in 
momentum space, the width in momentum space and the production time, respectively. For 
simplicity we assume that all the wave packets are emitted at the same instant and with the same 
width, = (vTj, ^j, 0", to)- The n boson states, normalized to unity, are given as 

I ai, ••• , ttn) = 51 n(««l«'^«) «n ••• (4) 

\^(n)i=l / 

Here a^'^^ denotes the set of all the permutations of the indexes {1, 2, n} and the subscript sized 
o-^ denotes the index that replaces the index « in a given permutation from a*-"-'. 

Solution for a New Type of Density Matrix There is one special density matrix, for which one 
can overcome the difficulty, related to the non- vanishing overlap of many hundreds of wave-packets, 
even in an explicit analytical manner. Namely, if one assumes, that we have a system, in which 
the emission probability of a boson is increased if there is an other emission in the vicinity: 



p„(ai,...,a„) = — ^J|Pi("i) n («fc|a<xfe) (5) 

■'^y'')i = l \^{n) fc = l / 

The coefficient of proportionality, M{n), can be determined from the normalization condition. 
The density matrix of eq. (H) describes a quantum-mechanical wave-packet system with induced 
emission, and the amount of the induced emission is controlled by the overlap of the n wave- 
packets yielding a weight in the range of [l,n!]. Although it is very difficult numerically to 
operate with such a wildly fluctuating weight, we were able to reduce the problem ^ to an already 
discovered "ring" - algebra of permanents for plane- wave outgoing states 0. 

For the sake of simplicity we assume a non-relativistic, non-expanding static source at rest in 
the frame where the calculations are performed. This may correspond to the central, mid-rapidity 
piece of the highly relativistic fireball created at RHIC. 



Pl(tt) = px{,i) Pp{,T^) 5{t - to) 
1 

1 

(27rmT)t 



P.(0 = 7:r^ exp(-eV(2i?')), 

Pp(7r) = ^ 3 exp(-7r7(2mr)), (6) 



and a Poisson multiplicity distribution p^"-' for the case when the Bose-Einstein effects are negli- 
gible: 

pi°^ = ^exp(-no), (7) 

This prescribes the multiplicity distribution of pion lasers in a very rare gas limit 0, and completes 
the specification of the model. The plane-wave model, to which the multi-particle wave-packet 
model was reduced in ref. [Q, can be further simplified to a set of recurrence relations with 
the help of the so-called "ring-algebra" discovered first by S. Pratt in ref. The probability of 
finding events with multiplicity n, as well as the single-particle and the two-particle momentum 
distribution in such events is given as 

Pn = UJni^UJk] , (8) 
\fc=0 / 



<^(ki) = (9) 

i=l 
n l-l 

iV(")(ki,k2) = ^[G^(l,l)G,_™(2,2) + G„(l,2)G,_a2,l)], (10) 

1=2 m=l 

where Un = Pn/Po and 



Gn{i,j) = n^/i„exp(-a„(k^2 + k]) +gnkikj). (11) 
Averaging over the multiplicity distribution p„ yields the inclusive spectra as 

oo 

G(l,2) = Y.Gn{l,2), (12) 

n=l 

oo 

iVi(ki) = ^p„iVf)(ki) = G(l,l), (13) 

n=l 

iV2(ki,k2) = G(1,1)G(2,2) + G(1,2)G(2,1). (14) 
With the help of the notation 

TTlT 

a^ = a^ + 2mT, Rl = + (15) 

a (Trp 

the recurrence relations that correspond to the solution of the ring-algebra [||, || are obtained 
for the case of the multi-particle wave-packet model. These correspond to the pion laser model of 
S. Pratt when a replacement R ^ Re and T ^ Te = crf./(2m) is performed. 
Let us introduce the following auxiliary quantities: 

7± = ^ (l + X ± Vl + 2x) X = Rial (16) 

The general analytical solution for the multiplicity distribution of the model is given through the 
generating function of the multiplicity distribution p„ 

oo / oo \ 

G{z) = Pn^"" = exp ^ Cniz- - 1) , (17) 

n=0 \n=l / 



where Cn is 



n 



7+ -7- 



;i8) 



together with the general analytic solution for the functions G'„(l,2): 



G„(l,2) = j„e 



n ri\2/n ri\2 

7jki-7jk2 + 7lk2-75ki 



TT 



1 7+ - 7- 
4 7!^-7- 



(19) 
(20) 



The detailed proof that the analytic solution to the muhi-particle wave-packet model is indeed 
given by the above equations is described in ref . . The representation of eq. (|l^ indicates that 
the quantities C^-s are the so called combinants |]TD|, [TT], |T2| of the probabihty distribution of pn 
and in our case their explicit form is known for any set of model parameters, as given by eqs. 
( [T^JT^JT5| ) . One can prove that the mean multiplicity (n) = J2'^=i'>^Pn = J^iZi^Ci. The large 

n behavior of nCn depends on the ratio of ^0/7+, since for large values of n, we always have 
(7-/7+) 2' << 1. The critical value of Uq is 



rir. 



3 

7| 



1 + X + Vl + 2x 



:2i) 



If Uq < Tic, one finds /zm„ 



cx,nCn = and (n) < 00, if no > ric one obtains linin-^cyonCn = 00 
and {n) = 00, finally, if = one finds linin^^nCn = 1 and (n) = 00. The divergence of the 
mean multiplicity (n) is related to condensation of the wave-packets to the wave-packet state with 
TT = 0, i.e. to the wave-packet state with zero mean momentum [Q, if no > n^.. 
The multiplicity distribution of eq. ([17|) is studied at greater length in ref. f^. 
Dense Bose-gas Limit: This wave-packet model exhibits a lasing behavior in the very dense 
Bose-gas limit, which corresponds to an optically coherent behavior, characterized by a vanishing 
enhancement of the two-particle intensity correlations at low momentum, C(ki,k2) = 1, a case 
which is described in greater details in ref. [Q. 

Rare gas limiting case. Large source sizes or large effective temperatures correspond to the 
X >> 1 limiting case, where the general analytical solution of the model, presented above, be- 
comes particularly simple and the exclusive and inclusive spectra and correlation functions can be 
obtained analytically to leading order in 1/a; << 1. From eq. ([TT]) one obtains that 



G„(l,2) 



3r. 



n 



(vra^)2 




i{n~l) 



[22] 



(23) 



where Ak = ki — k2. We can see from eq. (0), that the higher order corrections will contribute 
to the observables with reduced effective temperatures and reduced effective radii. 

The very rare gas limiting case corresponds to keeping only the leading n = 1 order terms in 
the above equations. The multiplicity distribution is a Poisson distribution with {n) = uq and 



no influence from stimulated emission. The momentum distribution is a Boltzmann distribution, 
and the exclusive and inclusive momentum distributions coincide. The leading order two-particle 
Bose-Einstein correlation function is a Gaussian correlation function with a constant intercept 
parameter of A = 1 and with a momentum - independent radius parameter of i?* = -Rg 0- 

The two-particle exclusive correlation functions can also be evaluated by applying a Gaussian 
approximation to the leading order corrections in the x >> 1 limiting case: 



n 



n[n 



iV(")(ki)iVr(k2; 



Ml 



1 + Akgxp 



(24) 



where K = 0.5(ki + k2), the side and outwards directions are introduced utilizing the spherical 
symmetry of the source as Ak^ = Ak - K(Ak ■ K)/(K ■ K) and Ak^ = K(Ak ■ K)/(K ■ K), 



similarly to refs. |13]. The mean-momentum dependent intercept and radius parameters are 



A 



K 




n + 2)Rl + — 



ex. 



T / 



(25) 
(26) 
(27) 



Thus the symmetrization results in a momentum - dependent intercept parameter Ak that starts 
from a Ak=o < 1 value at low momentum and increases with increasing momentum. Already in 
the first paper about the pion laser model, ref. |p, a reduction of the exact intercept parameter 
was observed and interpreted as the onset of a coherent behavior in the low momentum modes. 
First, a partially coherent system is created, characterized by Ak < 1, and if the density of pions 
is further increased, one finds a fully developed pion-laser with Ak = 0, see ref. for analytic 
considerations. 

Observe that the radius parameter at low mean momentum decreases while the radius param- 
eter at high mean momentum increases, as compared to Re- The radius parameter of the exclusive 
correlation function thus becomes momentum- dependent even for static sources! This effect is more 
pronounced for higher values of the fixed multiplicity n, in contrast to the momentum dependence 
of Ak that is independent of n. 

Last, but not least, a specific term appears in the two-particle exclusive correlation function, 
that contributes only to the out direction, which, in case of spherically symmetric sources, may 
be identified with the direction of the mean momentum [jl3|. This directional dependence is 
related only to the direction of the relative momentum as compared to the direction of the mean 
momentum, and does not violate the assumed spherical symmetry of the boson source. The effect 



vanishes both at very low or at very high values of the mean momentum K, according to eq. (|27|) . 

Highlights : Our new analytic result is that multi-boson correlations generate momentum- 
dependent radius and intercept parameters even for static sources, as well as induce a special 
directional dependence of the correlation function. The effective radius parameter of the two- 
particle correlation function is reduced for low values and enlargened for large values of the mean 
momentum in the rare gas limiting compared to the case when multi-particle symmetriza- 

tion effects are neglected. From the numerical studies, described in ref. 0, we find that for an 



extended, hot and rare gas of a few hundred pions, the reduction of the radius parameters at low 
momentum is found to be the most apparent effect. The directional dependence of the radius 
parameters and the enhancement of the radii at high momentum is characteristic for a small, 
cold pion gas with only a handful of particles in it. These results can be understood qualita- 
tively by an enhancement of the wave-packets in the low momentum modes, due to multi-particle 
Bose-Einstein symmetrization effects, as the system starts to approach the formation of a laser, 
characterized by the appearance of partial optical coherence in the low momentum modes. 
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